Electron dynamics in vacancy islands 
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The dynamics of Ag(lll) surface state electrons confined to nanoscale hexagonal and triangular 
vacancy islands are investigated using scanning tunneling spectroscopy. The lifetimes of quantised 
states with significant amplitude near the centers of the vacancies are weakly affected by the geome- 
try of the confining cavity. A model that includes the dependence of the lifetime on electron energy, 
vacancy size, step reflectivity and the phase coherence length describes the results well. For vacancy 
islands with areas in the range ~ 40-220 nm 2 lossy scattering is the dominant lifetime-limiting pro- 
cess. This result and a corrected analysis of published experimental data improve the consistency 
of experimental and calculated surface state lifetimes. 
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Spectacular real-space observations of electron con- 
finement have been realized by scanning tunneling mi- 
croscopy (STM) and scanning tunneling spectroscopy 
(STS) for the Shockley-type electronic surface states on 
the (111) surfaces of the noble metals. Electron con- 
finement to artificial (e.g., Refs. 1.2 3 4 5) and natural 
(e.g., Refs. 6,7,8 y) nanostructures has been studied in 
some detail. Moreover, several reports have addressed 
the lifetime of surface state electrons as deduced from 
tunnelling spectroscopy of large terraces, scattering pat- 
terns near steps and in confining structures ASiiSiiiiiS 
Currently, however, the influence of the confinement on 
electronic lifetimes is not clear. 

Using low-temperature STS of triangular and hexago- 
nal vacancy islands and model calculations, we investi- 
gate the effects of the geometry of the confining resonator 
and of lossy scattering at its boundary on the lifetimes 
of confined states. For the states investigated, which 
have an antinode at the center of these resonators, we 
find similar lifetimes independent of the geometry. Lossy 
scattering at the confining step edges turns out to be the 
lifetime-limiting process. 

The experiments were performed in ultrahigh vacuum 
(base pressure 1 x 10 -8 Pa) at 9 K using a custom-made 
microscope. Tungsten tips and Ag(lll) were prepared 
by argon ion bombardment and annealing. Hexagonal va- 
cancy islands were fabricated by briefly (2 — 3 s) exposing 
freshly prepared Ag(lll) surfaces to a low- flux argon ion 
beam. 13 The resulting nanostructures are shown in Fig. 
[I] Triangular vacancy islands were obtained through con- 
trolled tip-sample contacts. The differential conductiv- 
ity (dl/dV) was measured by superimposing a sinusoidal 
voltage modulation (3mV rms , 10 kHz) on the tunneling 
voltage and measuring the current response by a lock-in 
amplifier. 

Maps of dl/dV from a hexagonal vacancy island (Fig. 
Et) exhibit strongly voltage-dependent features, similar 
to those found from, e.g. hexagonal Ag islands. 6,8 For 
voltages below the Ag(lll) surface state band edge at 



FIG. 1: (a) Topographic image of Ag(lll) surface after low- 
flux Ar + bombardment, (b) Atomically resolved vacancy is- 
land which serves as the cavity for electron confinement. 



Eq « —67 mV the interior of the vacancy island is fea- 
tureless while at higher voltages increasing numbers of 
sixfold azimuthally symmetric rings occur. The sixfold 
modulation is partially distorted due to deviations of the 
vacancy island from a perfect hexagonal shape. dl/dV 
being related to the local density of states^ the inter- 
ference patterns are attributed to surface state standing 
waves confined by scattering at the edges of the island. 
Figure 03 shows similar data for a triangular vacancy is- 
land. Starting from a featureless interference pattern at 
— 60 mV, increasingly complex patterns occur at higher 
voltages. The number of antinodes increases with higher 
voltages as expected from a particle-in-a-box model. At 
— 20 mV, for instance, six antinodes of the local density 
of states are clearly visible. 

For a quantitative measure of the confinement we ac- 
quired dl/dV spectra above the vacancy island centers. 
A typical spectrum of a hexagonal island (Fig. |3Ji) is 
comprised of a series of peaks with the first appearing 
just above the lower band edge of the Ag(lll) surface 
state. The peak width increases with increasing energy. 
These peaks correspond to the energy levels of the con- 
fined Ag(lll) surface state electrons, with the broaden- 
ing into resonances as a result of single-particle scattering 
pro cesses pil^ many-body interactions^ and instrumen- 
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FIG. 3: (a) Experimental dl/dV data (dots) fitted by cal- 
culated spectrum (full line) based on a variation calculation 
as explained in the text, (b) Experimental dl/dV spectrum 
(dots) fitted by superposition of Lorentzians (full lines). To- 
pographic edge length of the island is 7.4 nm. 



FIG. 2: (a) Series of constant-current dl/dV maps of a hexag- 
onal vacancy island with edge length ^5.8 nm recorded at in- 
dicated sample voltages (white corresponds to large dl/dV). 
(b) Analogous data for a triangular vacancy with edge length 
22 nm. 



tal effects such as the finite modulation voltage and finite 
temperature. We have confirmed this interpretation us- 
ing the same approach as in Refs. IfillcSl using a variational 
embedding technique 17 to compute the electronic struc- 
ture of two-dimensional electrons with effective mass m* 
confined to hexagonal or triangular domains by an in- 
finite barrier. This gives a series of discrete levels at 
energies that are dependent upon the size and geometry 
of the confined region, which we broaden by including an 
empirically determined self-energy ImS « 0.2 (E — E-p), 
where Ep is the Fermi energy. Note that in compar- 
ing the experimental data with calculated spectra the 
true dimensions of the island must be determined. It is 
not a priori clear that the relevant dimension of the is- 
land coincides with its apparent topographic boundaries, 
for example due to the finite spatial extension of the tip 
the topographic dimensions are expected to be smaller 
than the actual ones. Consequently, we fit the experi- 
mental d//dV curves by calculated spectra varying the 
edge length as the only fit parameter. The full line in 
Fig. depicts the result of this fit procedure applied 
to the vacancy island with a topographic edge length of 
7.4 nm. The fit parameter, however, is given by 7.6 nm. 
Performing this fit procedure for several island sizes we 
find that the topographs of the vacancy islands tend to 
underestimate the actual sizes. 

We now discuss the lifetimes of the quantised electron 
states which are proportional to the inverse width of the 
associated peaks in the dl/dV spectra. In particular the 
half width at half maximum (HWHM) of the electron lev- 
els corresponds to an energy uncertainty (ImE) related 



to the lifetime r of the electronic state via 



r = ft/ (2 ImE) ~ [329 meVfs] /ImE. 



(1) 



In analysing the peak width, instrumental broadening 
due to voltage modulation and thermal effects have been 
taken into account. A fit using 5 Lorentzians plus an ad- 
ditional one at high energies to match the background is 
shown in Fig.^D. Taking the HWHM of the Lorentzians 
we arrive at the corresponding imaginary parts of the self 
energy, which is converted to a lifetime using fQ. Fig- 
ure^displays the electronic lifetimes versus their binding 
energies as extracted from various hexagonal and trian- 
gular vacancy islands, along with data from several other 
Ag(lll) surface state lifetime studies, both experimental 
and theoretical. While a general trend of increasing life- 
time as the binding energy approaches the Fermi level is 
obvious, there are significant differences in magnitudes, 
to be discussed next. 

We first note from the vacancy data of Fig. 0] that the 
geometric shape of the confining structure has little im- 
pact on lifetimes. Consequently, in modeling, we consider 
a circular vacancy island where the high symmetry per- 
mits a relatively simple description. In particular, one 
can show 18 that the confined levels occur at energies for 
which 
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(2) 



S is the radius of the vacancy, k = ^2m*(E — Eo)/h 2 is 
the magnitude of the surface electron wave vector, and 
R is the reflection coefficient describing scattering from 
the confining atomic step. Eq. [21 is similar to the round 
trip phase condition for surface states in the phase accu- 
mulation mo del ^2*211 where electrons undergoing multiple 
reflections at barriers with reflection coefficients R\ and 
R r separated by distance d are found at energies given 
by R\R Y exp 2\kd = 1. In the vacancy island the round 
trip corresponds to starting and finishing at the center 
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FIG. 4: Lifetime versus binding energy E — Ep found for 
Ag(lll) surface state electrons confined to vacancy hexagons 
(o), vacancy triangles (a), triangular corrals of Ag adatoms 
(V, adapted from Ref. |5|), circular and rectangular corrals 
of Mn adatoms (□, adapted from Ref. 4); lifetimes extracted 
from standing wave patterns near step edges (0, adapted from 
Refs. rrTIIT3) : photoemission data (x, from Ref. 

El); 

results of 

theory are depicted as a full and a dashed line (from Refs. l2ll 
andll3. respectively). 



and only involves one reflection — so just one reflection 
coefficient appears in (0 — and in the vacancy island 
the wave function must have an antinode at the center 
to be visible in STS, so that the round trip phase corre- 
sponds to an odd number of half wavelengths. This is the 
effect of the exp — iir/2 factor in ((2J. One final difference 
is that also takes into account many-body interac- 
tions through the factor containing the phase relax- 
ation length due to electron-electron (e-e) and electron- 
phonon (e-p) scattering. L$ is related by the group veloc- 
ity v g = hk/m* to the corresponding inelastic scattering 
lifetime tj: L<$> = v z r\. 

Writing R = expi (cj)R — iln Eq. is satisfied 
when 

<j> R -iln|i?| + 2kS + iS/L$ - tt/2 = 2irn (3) 

for n = 1, 2, . . ., which has solutions at complex energies 
E — ilmE. Assuming ImS <C (E — Eq), the real and 
imaginary parts of Eq.U3give the energy of the n'th level 
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and corresponding width 
Im E„. = — ^ 



In \R\ 
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(4) 



(5) 



where v g , \R\ and L$ are evaluated at the energy E n . 
Hence the lifetime of the surface state electrons confined 
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FIG. 5: Hexagon (o) and triangle (a) data from Fig. 0] rep- 
resented as Im E versus binding energy E — Ef . 



within the vacancy islands is given by 



(6) 



where the lifetime tr associated with lossy scattering at 
the atomic step is 



tr 
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(7) 



At 225 meV, Echenique and co-workers 12 find t\ « 62 fs 
in calculations that include the surface band struc- 
ture, and which treat e-e interactions within the GW 
approximation^ and use the full Eliashberg spectral 
function for the e-p interaction. 21 Biirgi et alA have mea- 
sured the reflection coefficient for surface state scatter- 
ing at ascending steps on Ag(lll) and report \R\ « 
0.23±0.07at 225 meV. Using Eq.0 (and m* =0.4m e ) n 
this gives tr ~ (9 ± 2) fs for an island with edge 7.5 nm, 
and from (0 a lifetime of r ~ (7.9 ± 1.5) fs. This is 
consistent with our measured lifetimes shown in Fig. 0] 
The calculated result is exemplary for the overall consis- 
tency with our measured lifetime data. We can therefore 
conclude that at this energy the electron states within 
the vacancy island are roughly 8 times more likely to de- 
cay via lossy scattering at the confining boundary than 
via inelastic e-e and e-p scattering. Repeating this anal- 
ysis at other energies, we find that the lossy boundary 
scattering is the dominant lifetime-limiting process for 
all islands that we have studied. They would need to be 
typically an order of magnitude larger in dimension for 
the inelastic e-e and e-p scattering to become dominant. 

Biirgi et ali find that \R\ decreases with increasing 
energy. Using their values we find that — \n\R\ varies 
approximately as [E — Eq) 1 / 2 . Since v g oc (E — Eq) 1 / 2 
and since LJ 1 is negligible with respect to In \R\/S in Eq. 
Odue to the dominance of the lossy-scattering mechanism 
the level width increases approximately linearly with en- 
ergy. Figure shows our results for the linewidths as a 
function of energy, which demonstrate this behavior. A 
similar variation in level width was previously inferred by 
Li et al& for surface state confinement to adatom islands, 
where a phenomenological self energy Im E = 0.2 (E—Eq) 
was found to match the observed level broadening. Biirgi 
have shown that the reflection coefficient for de- 
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scending steps, which surround adatom islands, varies in 
a similar manner to that for ascending steps. 
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In Fig. 21 various other lifetime data are shown. 
The Mn corral values obtained by Kliewer et aim (□) 
have been obtained by analysing experimentally de- 
termined spectral linewidths using the two-dimensional 
"black-dot" scattering model introduced by Heller and 
coworkers i The lifetime is obtained using Eq. Q where 
ImS is the electron self energy needed to bring calcu- 
lated spectra into agreement with measured ones. In 
effect, the lifetime is assumed to be due to e-e/e-p scat- 
tering, and hence to be identified as ti, as the scattering 
properties of the confining adatom array is incorporated 
via the multiple-scattering equations. However, the na- 
ture of the system is similar to the present one, with the 
total linewidth a sum of contributions due to lossy scat- 
tering and the e-e/e-p scattering. Hence the assumption 
that Mn adatoms act as "black-dot" scatterers means 
that the lifetime values that have been obtained should 
be viewed as lower limits to t\. The true reflectivity of 
the Mn adatom corrals is likely to be lower than that 
of the assumed "black-dot" scatterers, meaning that the 
lossy scattering is actually greater than modelled. Hence 
a smaller amount of the measured linewidth should be 
attributed to e-e/e-p scattering, increasing r\. 

The lifetimes deduced from the decay of standing wave 
pattern o 11 ! 12 at steps shown in Fig. 0] (0) are not compli- 
cated by lossy scattering as the asymptotic decay of the 
standing waves depends only upon r\. However, Biirgi 
et alA have misidentified the correct dependence of the 
standing wave patterns on the phase coherence length, 
so that their are actually twice t\? 2 The lifetimes de- 
termined by Braun and Rieder— shown in Fig. 0] (v) 
have been obtained from a detailed analysis of stand- 
ing wave patterns in triangular corrals constructed from 



Ag adatoms. In their analysis the Ag adatoms are not 
treated as "black-dot" scatterers, but as point scatter- 
ers whose scattering properties are determined by fitting, 
along with a phase coherence length which enters an at- 
tenuation factor exp— r/L& in the electron propagation 
between scattering events, and to and from the STM tip 
position. Thus lossy scattering effects are fully accounted 
for in this work - in so far as the two-dimensional point- 
scatterer model is valid - but again an incorrect rela- 
tionship is used between the phase coherence length L$ 
and the lifetime ti, leading to values that are twice t\. 
Correcting for this, the agreement between the lifetimes 
found by Braun and Rieder and theoretical values for 
Ag(lll), which are generally lower than those found in 
Ref. |3, is improved. 

In summary we have investigated the lifetimes of elec- 
trons confined to vacancy islands with areas « 40 — 
220 nm 2 on Ag(lll). We find that the geometry of the 
vacancy has only a weak influence on the lifetimes, which 
are dominated by lossy scattering at the steps which form 
the edges of the island. The lifetimes are well described 
by a theory that includes the dependence on the elec- 
tron energy, vacancy size, step reflectivity and the phase 
coherence length. This indicates that the crossover to 
lifetimes dominated by electron-electron and electron- 
phonon scattering will occur for islands with dimensions 
an order of magnitude greater than those studied here. 
This result and a corrected analysis of published experi- 
mental data lead to a more consistent picture of surface 
state lifetimes in experiments and calculation. 
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